Hidden XY structure of the bond-charge Hubbard model 
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The repulsive one-dimensional Hubbard model with bond-charge interaction (HBC) in the su- 
perconducting regime is mapped onto the spin-1/2 XY model with transverse field. We calculate 
correlations and phase boundaries, realizing an excellent agreement with numerical results. The 
critical line for the superconducting transition is shown to coincide with the analytical factorization 
line identifying the commensurate-incommensurate transition in the XY model. 
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The Hubbard Hamiltonian and its extensions are 
known to model several correlated quantum systems, 
ranging from high-T c superconductors to cold fcrmionic 
atoms trapped into optical lattices In particular, the 
HBC model describes the interaction between fermions 
located on bonds and on lattice sites 0, Q • This exten- 
sion is considered to be especially relevant to the field of 
high-T c superconductors [4|. In fact, it has recently been 
found [jjlfl that a superconducting phase takes place also 
for repulsive values of the on-site Coulomb interaction. 
The phase is characterized by incommensurate modula- 
tions in the charge structure factor. Its boundaries have 
been explored numerically, though their fundamental na- 
ture has not been understood yet. 

We find that the explanation of the above features re- 
sides into the underlying effective model, which for the 
superconducting phase turns out to be the anisotropic 
XY chain in a transverse field. Such model is known to 
be equivalent to free spinless fermions and it is remark- 
able how it can faithfully describe quantities of a strongly 
correlated system like the HBC chain. Indeed, the map- 
ping allows us to derive analytical expressions for both 
the critical line and correlations, reproducing with amaz- 
ing accuracy the numerical data. 

The model Hamiltonian for the HBC chain reads 



H = -^T / [1-X {n is + n i+lg )} (4 a c l+la + c\ +la c ia ) 

id 

U 
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where a =^, \. [a denoting the opposite of a), and the op- 
erator c ia creates a fermion at site i with spin er. More- 
over riic = cl^Ciu. The parameters U and X, expressed 
in units of the hopping amplitude, are the on-site and 
bond-charge Coulomb repulsion respectively. 

While the HBC model cannot be exactly solved for all 
A, there are two integrable point at X = and X = 1, 
for all values of U. The former is the well-known Hubbard 
model which is solvable by Bethe Ansatz. The integra- 
bility of the case X = 1 is due to the fact that the empty 



and the doubly occupied sites in this case are indistin- 
guishable, and the same holds for the t and 4- spins in the 
singly occupied sites, so that the model can be rephrased 
in terms of tight- binding spinless fermions in ID 0. In 
addition, the number of double occupancies turns out to 
be a conserved quantity. 

In the general case, Eq.flT]) can be fruitfully recasted 
passing to a slave boson representation. One can make 
the transformation |0) -> ej|0), cfjO) -> fj a \0) and 
cj^ctjjO) — > di\Q), where empty and doubly occupies sites 
are bosons, while the single occupations are fermions. 



The hard-core constraint eje^ + d\di + ^2 n<jfia = 1 
completes the identification. Then, the c-fermions are 



fLei + difia and n io 



flfir + djdi. The 



total number of particles is N — Nf + 2Nd- The filling 
factor is v = N/L, with < v < 2. Accordingly, we have 
v e + Vf + is c i = 1 and v = Vf + 2vd- After the substitution, 
the Hamiltonian becomes % = ^2 iG T~Lia, where 



7-Li a - -y/^/iCT + fi a fi+i,a (txd\ +1 di - e\ +1 ei 
- s xfhfi+i,a ( e i+i d i + di+ie.i) + H.c. 



(2) 



with tx = 1 — 2A, and sx = 1 — X. It can be rec- 
ognized that the first two terms describe the kinetic en- 
ergy of a single electron (hole) with spin a in a back- 
ground of empty (doubly occupied) sites, whereas the 
third term describes the transformation of two opposite 
spins into an empty and a doubly occupied site. Since 
the coefficient sx turns out to give the smallest contri- 
bution for X > 2/3, it is not surprising that the exact 
solution obtained assuming sx — (and arbitrary tx) 
Q shares in this regime many features of the ground 
state of the true model, obtained by numerical investiga- 
tion 0]. To some extent, these features hold within the 
range X > X c = 1/2, where X c is the value at which tx 
changes sign. This is true in particular as for the presence 
of phase coexistence of domains formed by only empty or 
doubly occipied sites, in which the single particles move. 
On the other hand, fixing sx = yields to a critical curve 
Ups = 4A for the stabilility of the phase separated re- 
gion, whereas the superconducting transition takes place 



2 



(only for sx ^ 0) at a value Use which is well below that 
line. Since the role of empty and doubly occupied sites, 
as well as the conservation of their number, appears to 
be the same as for sx = also in the superconduct- 
ing case, one can infer that it is just the motion of the 
single electrons and holes which determines the change 
Ups — > Use f° r s x 7^ 0. In this paper, we assume this 
point of view: treating the empty and doubly occupied 
states as the vacuum in which the single particles move. 

Let us go back to Eq.([2]) and consider what happens 
at sx 7^ 0. The SU (2) charge symmetry is broken down 
to U(l), which merely describes the conservation of the 
number of fermions. The large spin degeneracy is re- 
moved and it is like as if the fermionic dynamics is in- 
fluenced by the background imposed by the bosons and 
viceversa. This picture is correct as far as the spin and 
charge degrees of freedom are not separated. For X < 1, 
the pair creation term in Eq. d2J) induces short-ranged 
antiferromagnetic (AFM) correlations in both spin and 
pscudospin degrees of freedom. Since at half-filling the 
probabilities of having an empty and a doubly occupied 
sites are identical and coincide with 1/2, we can ap- 



proximate the term (ej, ^ — txd i+1 di) rj X. Thus, 
in this case the kinetic energy term in T-ii a becomes 
-X flf l+la + (X-l) flfl+i,* + Hc - where the term 
fiafi+i s al wa ys takes place due to the bosonic AFM cor- 
relations. Both the bosonic species are considered as a 
unique vacuum for the fermions /. 

Assuming the existence of AFM correlations also in the 
fermionic variables, we can drop the spin indices. The 
effect of f}f} + i is to open a gap at the Fermi level, hence 
reducing considerably the ground state (GS) energy. This 
mechanism is analogous to what happens in the case of 
the Peierls instability (in that case the gap is opened by 
the dimerization) where the bosons here play the role of 
the phonons that distort the lattice. So, we obtain a 
free-spinless fermion model TiS^ 



J2i=i H-i > where 
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It is instructive to notice that even in this form one can 
recover the exact solution of the case X = 1. Indeed 
a straghtforward diagonalization in Fourier space gives 
H = -2 J2k\- cos k + U/^flfk- The fermions fill the neg- 
ative energy states up to the Fermi point kf = i^Vf. The 
saturation occurs for U c = — 4cos(7Tf) for < v < 2. 

In the general case, W-f> can be easily shown to be 
equivalent to the following XY model in a transverse field 
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(4) 
at 



half filling. As usual we have applied the Jordan- Wigner 



transformation erf = 2fJ /j — I, o~f 

Ki-Ji with Ki = nLi (-°fc) = ex P[ i7r Efc=i n k\- 

AFM correlations in both bosonic and fermionic parti- 
cles are here assumed on the intuitive basic observation of 
the reduction of GS energy by means of the pair creation 
terms. A more rigorous approach would involve a self- 
consistent determination of the hopping coefficients in 
the quadratic model in which the spin labels are retained. 
Such approach allows to extend the analysis away from 
half filling and in magnetic field, and goes beyond the 
purpose of the present paper. We dedicate a forthcom- 
ing extended manuscript to a self-consistent approach. 

In what follows, we examine some important conse- 
quences that can be derived from the exact solution of 
the XY model, written in Eq.(|4]). 

As known, the Hamiltonian @ can be diagonalized: 

H = E Q + i J^keBZ^k (ptfik - I J 1 where the sum is 
performed in the Brillouin zone (BZ), and the dispersion 

relations are A^ = 1\j (cos k + h) 2 + 7 2 sin 2 k. Given 
the positiveness of A&, the GS energy Eqs is deter- 
mined by the vacuum of the Bogoliubov quasiparticles 
0k, giving E GS = E - jzJ2keBZ A k- B Y taking the 
thermodynamic limit L — > 00, we get an energy density 

We have compared the outcomes of our mapping with 
numerical calculations usin g th e density matrix renor- 
malization group (DMRG) [12(. In particular, we used 
extrapolations in 1/L of data collected by selecting 7 
finite-system sweeps and 1024-1152 states. Numerical 
and analytical results of the energy density at X = 0.8 
are displayed in table |U 
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-0.9717 
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0.173 


0.17588 


10/30 


0.33914 


2.5 


-1.3300 


-1.3287 


0.1063 


0.11488 


6/30 


0.20116 



Table I: Comparison between various quantities defined in 
the text computed either numerically (num) with DMRG or 
analytically by means of the equivalent XY model (th) for 
X = 0.8, both for periodic boundary conditions. The latter is 
treated directly in the thermodynamic limit, while the former 
are extrapolated to L — >• 00 from finite-size data. The char- 
acteristic wavenumber q is extracted from Fourier transforms 
at L = 30. 

An important feature of the XY chain is the presence 
of a factorization line h 2 +j 2 = 1, which corresponds to a 
commensurate-incommensurate (CIC) transition. In the 
HBC model this transition is mapped analytically into 

Use = 4V2A - 1. (5) 

Such transition was discovered numerically in Ref.Q 
and separates a incommensurate singlet superconducting 
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Figure 1: Comparison between the phase diagram of the HBC 
chain, calculated numerically in Ref. @] (symbols with dashed 
lines) and the phase diagram obtained from the mapping onto 
the XY model in transverse field (continuous lines). The up- 
per curves correspond to the spin gap transition where spin 
excitations become gapless, while the lower curves mark the 
transition into the ICSS phase where the charge compressibil- 
ity diverges. 



(ICSS) phase from a bond ordered wave (BOW) phase 
[6|. As seen in Fig|TJ the curve obtained with our map- 
ping describes rather accurately the numerical data of 
the transition. 

Along the factorization line the GS in the S = 1/2 
model is written as ®f =l \(f>), where |</>) = cos|| t) + 
sinf | |), with cos6» = [(1 - 7 )/(l + 7)] 1 / 2 = a. Here 
the local magnetization is 2vf — 1 = a = \J2X — 1. Ac- 
cordingly, the number of double occupations along the 
factorization line at half filling is v& = (1 — a)/4 = 
(1 — s/2X — l)/4. In the rest of the phase diagram, the 
transverse magnetization of the XY chain is given by 
= iY.keBz( h + cosfc)A~ 1 . The diverging charge 
compressibility of the ICSS phase is explainable simply 
by observing that adding two particles produce the con- 
version of an empty site onto a doubly occopied one, with- 
out changing the energy in the XY representation. 

In addition, the XY model in ID is known to undergo 
a quantum phase transition along line h = 1, belonging 
to the universality class of the classical Ising model in 
2D. This translates directly into the line U = AX in the 
phase diagram of the HBC model (see FigfT]). The latter 
coincides with the critical line of stability of PS (Ups) in 
the integrable case sx — 0, and is close to the numerical 
critical line between the spin density wave (SDW) and 
the BOW phase in Fig. 1, at least for X close to 1. 
Moreover, the line 7 = 1, which is known to describe 
the Ising model in a transverse field, here corresponds to 
the case X = 1/2. While it is questionable whether the 
assumptions that have originated our approximations for 
the ICSS phase are still valid in the above limiting cases, 
one can recognize that instead at the very crucial critical 



U= 2.5, X= 0.8 




Figure 2: Analysis of the various contributions [see Eq.Q] to 
the static structure factor (Fourier transform) of the density 
correlation function for the HBC model [L = 30) with the pa- 
rameters reported in the legend. The vertical line corresponds 
to tfi/ir (see tabled]). 

point X = 1/2 and U = 0, our system described in Eq.(|TJ) 
is mapped into nothing but the Ising model. 

From the seminal paper of Barouch and McCoy on 
the statistical mechanics of the XY model [ll|, it is 
known that the oscillation wavenumber of the correla- 
tor Pxx{R) = (cfCi+p) m t ne incommensurate region 
h 2 + 7 2 < 1 is 

Jl - 7 2 - h? I2X - 1 

tanw= ; = a 5- - 1, (6) 

h V (C//4) 2 

with a period Rq = 2ir/ip. 

A first striking observation is the fact that the correla- 
tions of the total density exhibit a peak very close to the 
characteristic wavenumber ip in Eq.©: in the last two 
columns of table U we report the wavevector q at which 
the total density structure factor has a peak (see an ex- 
ample in fig. [5]) and the corresponding value of "0. For 
X = 0.9 and U = 3 with L — 32 the peak is located at 
q/n = 0.1875 while ip/ir = 0.18342. 

The appearance of the peak at wavenumber Q in the 
Fourier transform of a correlation function that decays 
as cos(QR) cxp(— R/£)/R a is related also to the expo- 
nent a: the smaller is a the sharper is the peak. In 
particular for a = 2 which is the case for the correla- 
tion p zz {R) = (ofof +i? ) of the XY model, the peak it 
not visibile at all, despite the fact that the oscillations 
actually have characteristic wavenumber 2if>. Hence, it is 
worth to inspect in more detail the origin of the peaks 
observed numerically. Since the local density operator m 
in terms of single and double occupancies n S i and is 
given by rii = n S i + 2n,a , the correlation function of the 
total density decomposes in the following parts 

(mn i+R ) =(n si n 3i+R ) + 4(n di n dt+ p) 

+ 2{n dl n sl+B ) + 2(n si n di+R ). (7) 
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Figure 3: Comparison of the connected real-space correla- 
tion functions N S (R) (at half-filling) and p(R) for the HBC 
and XY model, respectively (see text for definitions based 
on singly occupied sites operators). The parameters of the 
two models are related by mapping as 7 = (1 — X)/X and 
h = U/4X. All the DMRG calculations for the HBC model 
and the analytical experessions of the curves for the XY model 
refer to L — 50. From top to bottom the data have been offset 
by +0.15, +0.10 and +0.05 for the sake of clarity. 

It turns out that the peak in the static structure factor 
is not due to the first term, but it is instead provided by 
{ndiUdi+Rj , as shown in Fig(3]for the test case X = 0.8 
and U = 2.5, although we obtained the same qualitative 
picture at U = 1. 

According to our mapping, we can compare directly 
the connected correlator N S (R) = (n S in S i +R ) — (n S i) 2 
in the HBC model with the density correlation function 
P(R) = (fi fifi +R fi+R) - (fi fi) 2 for the spinless fermions 
with Hamiltonian The calculation of the latter is 
omitted here since it is quite lengthy, though it simply 
involves a standard application the Wick theorem. The 
fully fermionic correlator N s (R) and the spinless fermions 
correlator p(R) are compared in Fig(3]for various choices 
of the parameters U and X in the ICSS phase of our 
starting system (i.e. the incommensurate one in the XY 
model); the agreement in real space is generally very 
good. Such behaviour of N S (R) is not obvious a priori 
in the HBC model and we interpret it as a remarkable 
nontrivial prediction of our mapping. 

In summary, we have studied the Hubbard model with 
bond-charge interaction in the superconducting regime, 
unveiling its underlying XY structure. We have shown 
that at half filling the numerical critical line for super- 
conductivity coincides with remarcable accuracy to the 
analytical factorization curve that marks the CIC tran- 
sition of the anisotropc XY model in a transverse field. 
Exploting the mapping for the calculation of correlations 
in the effective model has allowed us to predict rather 
accurately the peak in the charge structure factor of the 
original model. The results confirm a posteriori the cru- 
cial role of short range AF correlations and spin degrees 



of freedom as to the onset of superconductivity. The ul- 
timate presence of the latter is however to be ascribed to 
the interplay of the spin with the charge degrees of free- 
dom, the superconducting properties being absent from 
the incommensurate phase of the free fermions model. 

Based on the success of the present mapping, a number 
of further result are now in order. First, since the one- 
dimensionality of the model is not crucial to the map- 
ping, the latter should hold in higher dimension as well. 
In 2D, the numerical investigation of the XY model has 
been largely explored in the literature: this could pro- 
vide useful hints on the type of phase diagram which 
characterizes the 2D HBC model. Moreover, it would be 
interesting to understand the implications on the HBC 
model of a non-vanishing string order parameter which 
is peculiar of the XY model in transverse field. Finally, 
we expect that a similar mapping should hold also in the 
strongly repulsive regime U — > 00, since in that case no 
doubly occupied sites occur, and it is still quite natural 
to assume short range AFM order of single particles. 
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